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Precede de codage d*un maillage source par minimisation de 
volume, et applications correspondantes. 

L domaine de V invention et applications 
hi domaine de V invention 

5 Le domaine de T invention est celui du codage de stmctures de donn6es 

g6om6triques, ou maillages, notaniment de grande taille. Plus pr6cis6ment, 
rinvention conceme la representation et le codage d'objets ou de scenes en trois 
dimensions. Plus pr6cis6ment encore, T invention concerne une technique 
d' approximation d'un maillage source en trois dimensions. 

10 Un maillage est classiquement defini par un ensemble de sonmiets et de 

faces (ou un ensemble de sommets, d* aretes et d' orientations) d^finissant une 
topologie. De tels maillages sont par exemple utilises en graphisme sur ordinateur, 
pour mod61iser des objets en trois dimensions avec une complexity g^ometrique 
limit6e. 

15 L' approximation d'un maillage M consiste ^ trouver un maillage M' dont la 

complexite geometrique est moindre que celle du maillage M, et qui approche au 
mieux la g6om6trie de M, au sens d'un certain crit&re d6£ini Ji Tavance. 
i.2 exemples d* applications 

L* invention trouve des applications dans tons les domaines ou il est 
20 souhaitable de i^duire le nombre d' informations necessaires pour repr6senter et/ou 

manipuler efficacement un objet en trois dimensions ou un ensemble d'objets, par 
exemple pour r analyser, ou stocker et/ou le transmettre. 

A titre indicatif, 1' invention pent notamment s'appliquer au domaine de : 

la realite virtuelle (visites ou boutiques virtuelles, loisirs, 
25 telemanipulation, etc.). Dans ce type d' application, 

r approximation de maillages permet de reduire le coflt de rendu de 
scenes complexes, en particulier en d6finissant la notion 
d'echelonnabilite sur les maillages (fonction du point de vue, des 
capacit6s graphiques, du taux de rafraichissement souhait6, ...). 
30 Dans le cas de la r6alit6 virtuelle distribu6e ou partag6e, cela permet 
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6galement d' adapter la complexite d'une scene aux capacit6s de 
rendu et de stockage des differents terminaux, ainsi qu'aux debits 
des r6seaux ; 

la simulation scientifique (616nients finis, CAO, etc...)- La 

5 reduction de la complexite g^ometrique des modeles permet une 

acceleration des temps ae calcui, une prise de decisions plus^ 
rapide, notamment lors de la conception en CAO, et T Elimination 
des informations redondantes dans une base de donnees 3D ; 
la modfilisation (scanner 3D (reconstruction de surfaces k partir de 

10 points non organises), scanners volumiques, reconstruction de 

surfaces h partir de photos stereoscopiques ou de sequences video, 
de modules numeriques de terrains (imagerie satellite ou radar), 
etc.). Un module numerique de terrain permet ainsi Tobtention 
d'un maillage repr^sentant la topologie d'une region. Un tel 

15 maillage est obtenu par Techantillonnage regulier d'une image 

stockant T information d' altitude en chaque point. II en resulte une 
quantity de donnees importante, comprenant des informations 
inutiles pour la simulation scientifique, ou trop coflteuse pour le 
rendu (dans le cas de simulateurs). L' approximation de maillages 

20 reduit la quantite de donnees, tout en garantissant une bonne 

fidelite geometrique aux donnees initiales et la conservation de la 

topologie. [ 

2. art anterieur 

2.1 les families d'algorithmes 
25 Plusieurs techniques d* approximation de maillages sont deja connues. Les 

plus repandues peuvent Stre classees en trois grandes families d'algorithmes selon 
qu*elles fonctionnent par : 
decimation ; 

re-&hantillonnage sous-critique ; 
30 - subdivision adaptative. 
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2. LI decimation 

La decimation consiste a retirer de manifere iterative des sommets et/ou des 
faces d'un maillage. Cette operation est appelde operation 616mentaire de 
5 simplification. Les methodes mettant en oeuvre ce principe de decimation peuvent 

egalement optimiser les positions des sommets apres simplificauon, cene demie^ 
6tant choisie de maniere h, preserver au mieux la topologie du maillage. 

2.7.2 re-echantillonna^e sous-critique 

Le re-echantillonnage consiste k 6chaiitillonner un modfele original, soit en 
10 prenant des points aleatoirement sur sa surface et en retriangulant ensuite soit en 

definissant une grille tridimensionnelle et en agglom^rant les sonmiets dans chaque 
boite 61ementaire de la grille. Le mod^e ainsi g€n€r€ est simplifi6, et doit 
approximer au mieux les donntes initiales. Cette technique est rapide, mais ne 
conserve pas la topologie ni les caracteristiques visuellement importantes des 
15 maillages. 

2.7.3 subdivision adaptative 

La subdivision adaptative conmience avec un module comportant une 
geometric trfes simple, que Ton subdivise recursivement ensuite, en ajoutant h, 
chaque iteration un detail dans les regions ou Terreur d' approximation est 
20 maximale. 

2.2 exemples de techniques connues 

De fa9on k permettre une appfOXlmation d'un maillage avec une qualileTte" 
reconstruction satisfaisante, il est necessaire de combiner une decimation et une 
optimisation des positions des sonunets conserves. En d'autres termes, Tobjectif 
25 de base d'un proced6 de codage d'un maillage source etant de maximiser la qualite 

de r approximation pour une complexity geom6trique donn^e, celui-ci doit 
notamment presenter les proprietes suivantes : 
decimation ; 

preservation de la topologie ; 
30 - optimisation des positions, suivant un critdre d'erreur pr6defini. 
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^.2.1 renudUa^e 

Une premi&re m^thode , appel6e "remaillage" rcpondant h. ces critferes a et6 
pr6sent6e par Greg TURK, dans le document "Re-tiling polygonal surfaces" 
5 (SIGGRAPH 92 Conference Proceedings, pages 55-64, 1992). Elle fonctionne 

par fcchaiitillonnage, dtU;iiiiaUoii el upliimi>ailuu act> pomUuuh. EIIl pumiLi Ju 
param6trer le nombre de niveaux de resolution souhait6s et le nombre final de 
sommets pour chacun d'eux. 

L'algofithme debute par une distribution al6atoire du nombre de sommets 
10 param6tr6s sur la surface du modele, puis optimise leur position en les liant par 

des forces de repulsion (fonctions de la courbure locale). II triangule ensuite les 
polygenes ainsi formes, en incluant les sommets du maiUage original. Enfin, on 
dicime le module un sommet aprfes I'autre (si cette operation conserve la topologie 
locale) puis on retriangule le trou resultant. 
15 On obtient ainsi un module raffine dans les regions de forte courbure et 

d^cime sur les regions planes. 

Un inconvenient de cette metiiode est qa'elle pr6sente une partie al6atoire. 
En consequence, deux decimations successives ne donnent pas le mfeme resultat. 
Par ailleurs, cette technique a tendance h. lisser les discontinuites dans la plupart 

20 des situations. 

De plus, la fid61it6 aux donn6es initiales n'est pas exprimee directement au 

cuuis de ruptlniisatiou des positions. 

2.2.2 codage pro ^ressif de maillages 

Une autre technique, appelee "codage progressif de maillages", a 6t6 
25 d6velopp6e par Hugues HOPPE, dans le document "Progressive meshes" 

(SIGGRAPH 96 Conference Proceedings, pages 99-108, 1996). EUe repose sur 
la decimation et 1' optimisation des points. Cette technique est decrite plus en d6tail 
en annexe 1, en relation avec les figures 1 a 4, de fa$on h ne pas surcharger la 
presente discussion, n est clair cependant que cette annexe 1, ainsi que 1' annexe 2, 
30 font partie integrante de la description. 



Cette technique de codage progressif de maillages pr6sente differents 
inconvdnients. Notamment, elle ne conserve pas naturellement les discontinuit6s 
ou singularit6s geom6triques, ce qui suppose une gestion de cas particuliers, tels 
que les coins, les aretes vives finissantes et les arfetes vives r6guliSres, pour 6viter 
des ruptures de topologie. Cela entralne des calculs complexes, difficiles k mettre 
en oeuvre, et suppose des parametrages peu aises ^qu esi-ce qu unc uTei^ 



Vive 



Comment ponderer le terme de ressort pour 6viter de lisser les surfaces ? ...)• 

Par ailleurs, elle n'est pas efficace pour caract6riser les diff6rences 
perceptuelles entre deux maillages. 
10 nhiectifa de I 'invention 

L' invention a notamment pour objectif de pallier les differents 

inconv6nients de la technique anterieure. 

Plus pr6cis6ment, un objectif de 1' invention est de foumir un proc6d6 de 
codage d'un maillage source en trois dimensions, presentant un meilleur rapport 
15 quality de 1' approximation/complexite geometrique que les techniques connues. 

Notamment, 1' invention a pour objectif de foumir un tel proc6d6, 
lespectant plus fidMement les caracteristiques importantes k conserver (singularit6s 
ou discontinuit^s). 

Un autre objectif de I'invention est de foumir un tel proc6d6, ne n6cessitant 
20 pas une recherche et une gestion des cas particuliers (tels que les arStes et les 

coins), ni la mise en oeuvre de parametres specifiques (par exemple critferes de 

reconnaissance d'une arSte, seuUs, ...) qu'il est fl^eessaire Ue d^riiiii poui cliaquc 

maillage. 

L'invention a encore pour objectif de foumir un tel proc6de, qui permet 
25 vme reconstraction successive du maillage, celui-ci etant reconnaissable, sous une 

representation grossifere, trfes rapidement. 

L'invention a 6galement pour objectif de foumir un tel proc6d6 permettant 
de repondre aux fonctionnalit6s fixees dans le projet de standardisation ISO- 
MPEG4, a savoir : 

30 - la compression de maillages, qui doit r6duire fortement la quantity 



de donn^es d6crivant la geom6trie, avec des pertes visuelles 
minimales ; 

la representation progressive d*un maillage, permettant de definir 
un maillage sous la forme d'un maillage de base h faible resolution, 
5 et d*une sequence de raffmements ; 

r adaptation k une gamme de debus imporiiuUc, la 4Uttlilc cLmit 
acceptable, mais pour de tres faibles debits. 
4. caracteristiques principa l^x de V invention 
4,1 utilisation d'wie met rique volumique 
10 Ces objectifs, ainsi que d'autres qui apparaitront plus clairement par la 

suite, sent atteints selon F invention k I'aide d'un procede de codage d*un maillage 
source (M) representatif d'un objet en trois dimensions, d61ivrant un maillage 
simplifi6 (M') correspondant audit maillage source (M), lesdits maillages 6tant 
defmis par un ensemble de sommets, d'aretes et/ou de faces, Ce procede met en 
1 5 oeuvre une etape de minimisation d'un volume compris entre ledit maillage source 

(M) et ledit maillage simplifie (M'). 

Ainsi, r invention conceme une technique d' approximation de maillages 
utilisant une metrique basee volume, et non conmae cela est classiquement le cas, 
une metrique considerant la distance entre un sommet et une surface (une telle 
20 distance n'etant d' ailleurs pas unique). 

Cette approche nouvelle prend en compte natureUement les caracteristiques 

importantes A conserver (singularitds) sans qu'il suit ncccssairc de rccourir a la 

detection de cas particuliers. Aucun param6trage particulier n'est requis, avant le 
traitement d'un maillage : on minimise simplement le volume entre les deux 

25 maillages. 

4,2 caracteristiques avantageuses de Vinvention 

L' invention pent notamment €tre mise en oeuvre k Taide d'un algorithme 
organise en deux parties : la decimation et roptimisation de positions. 

La decimation permet d'obtenir un maillage geometriquement simplifi6, k 
30 partir du maillage original, tout en conservant la topologie et une bonne 
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ressemblance avec ce maillage original. 

L' optimisation des positions des sommets permet Tobtention de la 
meillexire approximation au sens de la metrique bas6e volume de 1' invention. 

En combinant decimation et optimisation, on pent g6n6rer soit un maillage 
5 sous fomxe progressive (en utUisant altemativement la decimation d'un sommet et 

une optimisation locale), soit des niveaux de resomuon aisimcis (,en uiilisimi 
successivement la decimation d'un ensemble de sommets et une optimisation 
globale). 

De fa^on avantageuse, chacun desdits maillages etant defini par la position 
10 de chacun de ses sommets et des aretes les reliant, ladite etape de minimisation 

assure la determination des positions des sommets (Xj, X2,-.. X^) dudit maillage 
simpUfie (M') minimisant le volume V(M, M') entre ledit maillage source et ledit 
maillage simplifie. 

Pref6rentiellement, ladite etape de minimisation met en oeuvre un 
15 processus it^ratif optimisant progressivement les positions des sonmiets dudit 

maillage simplifie (M'). 

Ledit processus iteratif pent notamment Stre interrompu lorsqu*au moins 

un des criteres d' arr§t suivants est atteint : 

un nombre (N) maximum d' iterations ; 
20 - une difference entre deux vecteurs deplacement successifs des 

positions des sommets inferieure k un seuil predetermine (e). 

Selon un mode de realisation avantageux, ladite etape de minimisation met 
en oeuvre une methode de gradient adaptatif. Cette methode, ainsi qu'on le verra 
par la suite, permet de mettre en oeuvre I'invention avec une complexite de calcul 
25 reduite. En effet, il n'est alors pas necessaire d'effectuer un calcul explicite du 

volume. 

Dans ce cas, ladite methode de gradient adaptatif peut avantageusement 
reposer sur les operations smvantes : 

selection d'un vecteur X de R (n > 1) dudit maillage simplifie. 



et calcul du gradient VE (X^) en Xp de la fonction h minimiser 
E = d(M,M'); 

determination de la position Xp* de Xp dudit maillage, selon la 
relation definie a I'iteration k+1 par : 



l|VE(X/)|| 



k variant de 0 a n-1 (avec n < N) et etant le pas de ladite 
relation ; 

N represente le nombre d' iterations utiles pour atteindre la 
convergence (cette demiere etant atteinte lorsque la difference entre 
10 deux determinations successives est inferieure ^ un seuil (e) ou 

lorsqu'un nombre maximum d' iterations est atteint). 
On pent noter que Xp appartient a si on effectue une optimisation 

locale ou k R^n si on effectue une optimisation globale. En d'autres termes, on 
optimise les positions d'un sommet a la fois dans le cas d'une optimisation locale 
15 et de tous les sommets simultanement (sous la forme d'un unique vecteur de R ^n) 
dans le cas d'une optimisation globale. 
De fagon prtferentielle, ledit pas jy est variable, et il varie en fonction de 

r oscillation de deux deplacements successifs des sommets et/ou en fonction des 
variations de Tenergie. II pent notamment etre maintenu entre deux homes et 

20 Y . 

' max 

Selon un mode de realisation avantageux de I'invention, on determine, h, 
chaque iteration, une variation elementaire dudit volume correspondant a un champ 
de vecteurs . La surface 6tant parametric par u et v, de fa9on qu'un vecteur 



SSt s'exprime sous la forme (u, v), ladite variation ^lementaire est assimilfe 
au paralleMpip&de g6n6r6 par revolution de T^lement de surface dudv dans la 
direction SST (u, v). 

Cette variation 616nientaire de volume entre M et M' pent 
""5 avantangeusement s'ecnre : 

SdCHM* ,5K* ) = - J J ii(u,v)n(u,v)8Sl' (u,v)da(u,v) 

u,v 

avec: ti(u,v) = ti(M,M ,n(u,v)) = 1 si la normale unitaire n(u,v) est 

orient^e vers le volume int6rieur, et -1 sinon. 
Selon un mode de realisation pr6f6rentiel de T invention, ledit maillage 
10 simplifi6 est param6tr6 k Taide d*un modfele d' Elements finis. Lesdits Elements 

fmis sont avantageusement obtenus k Faide d'un interpolatevir affijie. 
Dans ce cas, ledit maillage simplifi6 peut s'&rire : 

N 

M(u,v)= j:XiA.i(u,v) 

avec X..(u,v) une fonction de forme adaptee h un modele d'616ments finis 

15 triangulaires, 

et Xj sonmiet dudit maillage, d6fini dans R^, 
et on d^finit une fonction de forme sur le maillage par ie biais ae cooraonnees 
barycentriques. 

En utilisant I'expression pr6c6dente, on obtient : 



20 



N 

5d(M,M ,5Xi 5X J = - J J Ti(u,v)ii(u,v) ^ 6XiA,i(u,v)da(u,v) 

u,v *=1 

N 

= -LJJ Ti(u,v)n(u,v)5XiXi(u,v)da(u,v) 
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De fa9on avantageuse, la position d'un sommet Xj & la k*^™* iteration 

s'ecrit : 



ax, 



avec la derivte partielle de la distance pour un sonunet Xj : 

— — = - J J Ti(u,v)n(u,v)Xj(u,v)da(u,v). 

^Xj u,v eSupp(X^ 

Avantageusement, le precede de T invention met en oeuvre un codage 
progressif dudit maillage simplifie, par decimation et optimisation locale, Ainsi, 
ledit maillage simplifie est, de fagon avantageuse, represente par un maillage de 
base et une sequence de raffinements dudit maillage de base, 

Dans ce cas, le procede comprend preferentiellement une etape 
d' initialisation du vecteur Xp, comprenant une etape de denombrement des aretes 

vives autour de deux sommets k fiisionner, apres ladite decimation, et : 

si les nombres d'arStes vives sont les mSmes autour des deux 
sommets, on place le sommet resultant de la fusion au milieu du 
segment reliant lesdits sommets ; 

si les nombres d'argtes vives sont differents, on place le sonmiet 
resultant de la tusion sur le sonunet pfisentant le plus grand 



nombre d* aretes vives. 
De fa9on avantageuse, le calcul du gradient a la k^^"^« iteration comprend 
les 6tapes suivantes : 

discretisation de Texpression de la derivee partielle de la distance 
pour chaque sonamet Xp sous la forme : 

ddd " " 

= - Z Z Z Ti(i„jjn(i^,j,)>.i(i„j,)da(i„j,) 

TeSi=lj=l 



11 

avec : - S T ensemble des triangles voisins du sonunet . 

N-^ le nombre de points 6chantillonnes dans la 

direction de u ; 

le nombre de points ^chantillonnes dans la 

Hir^ ction d e v ; 



calcul de T orientation des surfaces, par identification de la plus 
proche intersection avec le maillage source M. 
Selon un aspect pr6ferentiel de I'invention, rorientation relative des 
surfaces audit maillage source et dudit maillage simplifi6 est calculee selon 
10 Tequation : 

Ti(u,v) = -<nM.,nM>.<nM., X^X^ > 

avec : Xj le point 6chantillonne sur M' ; 

Xj^ le point d' intersection de la droite passant par Xj et de direction n 

avec le maillage source M ; 
15 a M la nomiale au maillage source M au point Xj^ ; 

n,^. la noraiale au maillage source M' au point Xj ; 



( XjX j4 repr^sente un vecteur). 

Selon un autre aspect avantageux de I'invention, le proc6de comprend une 

^tapf> Hp. lin iitafmn Hp la degradation due h une transfomiation elementaire, mettant 

20 en oeuvre une queue de priority sur les transformations elementaires. 

Pref^rentiellement, ladite etape de limitation de la degradation due a une 
transformation €lementaire comprend les Stapes de : 

calcul d'un coiit pom chaque transformation el6mentaire possible ; 
realisation de la transformation 61ementaire de moindre cout ; 
25 - recalcul des couts des transformations dlementaires modifiees par la 

transformation 616mentaire pr&ddente ; 

ajout des nouvelles transformations 616mentaires cr66es et calcul 
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des coflts correspondants. 
Le coflt d'une transformation elementaixe (T^) pent par exemple s'exprimer 



par : 



C(T,(X„X.)) = max d,(VM,F(X[)) 
n^rr ' T] trnn*^^^"^^^^^" fn^innn ant deux sommets X; et Xj dudit maillage 



simplifieM' ; 

X.^ le sommet dudit maillage simplifie M' resultant de ladite 
transformation ; 

F(X j^) les faces dudit maillage simplifie M' avoisinant le sommet X.^ aprfes 

10 ladite transformation ; 

Vjvi ensemble des sommets dudit maillage source M appartenant aux 

faces ay ant ete intersectees lors du calcul de T orientation des 
surfaces pendant ladite minimisation. 
Le procede de T invention pent s'appliquer k de tres nombreux domaines 
15 techniques, et en particulier a au moins un des domaines appartenant au groupe 

comprenant : 

la realite virtueUe ; 

la simulation scientifique ; 

la modelisation. 

20 5. liste des figures 

D'autres caracteristiques et avantages de T invention apparaltront plus 
clairement a la lecture de la description suivante d'un mode de realisation 
preferentiel de I'invention, donnee a titre de simples exemples illustratifs et non 
limitatifs, et des dessins annexes parmi lesquels : 
25 - les figures 1 a 4 sont discutees en annexe 1, presentant la technique 

connue du codage progressif de maillages, et illustrent 
respectivement : 

figure 1 (egalement commentee par la suite) : 
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transformation el6mentaire (fusion d' arete) ; 
figure 2 : solutions initiales ; 

figure 3 : mise & jour des transformations dans le voisinage 

affect6 par la demifere modification ; 
5 - figure 4 : ajout des transformations nouvellement cr66es ; 

la figure 5 illuslie le pilncipe dfc base de T invention, cu^sisl^^ll ii~ 
minimiser le volume entre deux surfaces ; 

les figures 6 i 1 1 sont discut^es en annexe 2, decrivant les aspects 
mathematiques de I'invention, et illustrent respectivement : 
10 - figure 6 : variation 61ementaire du maillage ; 

figure 7 : variation elementaire du maillage dans la direction 
de la normale ; 

figure 8 : systeme de coordonnees sur un triangle ; 

figure 9 : f onction de forme ; 
15 - figures 10 et 11 : exemples d' Evolutions de surface sur 

deux courbes en deux dimensions ; 
la figure 12 illustre la mise k jour des adjacences de face apr&s 
fusion d'une arete adjacente a deux faces ; 

la figure 13 illustre la mise a jour des adjacences de face apres 
20 fusion d'une arete adjacente k ime face ; 

la figure 14 illustre le denombrage des aretes vives autour de 
sommcts du maillag e-; ■ 



la figure 15 presente le principe de T initialisation de la position du 
sommet a optimiser, dans differentes situations ; 
25 - la figure 16 illustre Techantillonnage sur un triangle pour 

r integration numerique ; 

la figure 17 explique le principe de Torientation relative des 
surfaces, de fafon k Eviter la formation de plis sur la surface 
losrque la position initiale est trop 61oign6e de Toptimum ; 
30 - la figure 18 illustre le coflt d*une transformation 616mentaire ; 
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la figure 19 est un synoptique de T approximation de maillages 
selon rinvention ; 

la figure 20 est un synoptique de 1' optimisation d'un somnaet ; 

les figures 21 k 24 illustrent le comportement du procede de 
5 rinvention, sur un parallelepipede rectangle, apres decimation sur 

un coin (figure 21), une arete vive reguliSre (figures '2:i et 23) et un 
plan (figure 24) ; 

la figure 25 illustre les diffSrentes (12) iterations effectuees sur un 
exemple simple, selon rinvention. 
10 6. principe general d'une mise en oeuvre de V invention 

6.1 structure 

Ainsi que discute precedemment, rinvention concerne une m6thode 
d' approximation de maillage mettant en oeuvre une metrique basee volume. Un 
des avantages de cette metrique est qu'elle prend naturellement en compte les 
15 caract6ristiques importantes h conserver (les singularites), sans devoir recourir a la 

detection de cas particuliers. 

D*une fafon g6nerale, T approximation de maillages genere un mod&le 
polygonal geometriquement plus simple que le modele original. Le nombre de 
polygones necessaire & la representation est ainsi reduit, tout en conservant la 
20 meilleure approximation, au sens perceptuel, du modele original. 

Le probleme correspondant pent gtre formule par une approche 
variationnelle permettant d'approximer un maillage M par irn maiiiage jvi' {on Ivl) 

comprenant un nombre reduit de triangles. Elle autorise le deplacement des 
sommets en optimisant leurs positions au sens d'un ou plusieurs criteres definis 
25 par le biais d'une fonctionnelle d'energie k minimiser. 

L' approximation d'un maillage genere un modfele geometriquement plus 
simple que le modele original. Le nombre de polygones necessaires a la 
representation est ainsi reduit tout en conservant la meilleure approximation (au 
sens perceptuel) du modele original. La qualite de T approximation definit le niveau 
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de ressemblance perceptuelle pour une complexity gtomdtrique fixee. L^invention 
utilise une m6trique bas6e volume entre "le maillage approximant" et le "maillage 
approxim6". 

Le principe de ralgorithme est scindd en deux parties : la decimation et 
5 r optimisation des positions. 

L' operation eiemeniaire de aecimauon chuisic csl la fuMuii d'aieie, qui = 
fusionne deux sommets, de fagon k supprimer un sommet et deux faces du 
maillage. L'algorithme d' optimisation choisi est un algorithme it6ratif de gradient k 
pas adaptatif afin de minimiser iterativement le volume compris entre le modele et 
10 le maillage original. 

En combinant ddcimation et optimisation, on peut g^n^rer soit un maillage 
sous forme progressive, soit des niveaux de resolution distincts. 
6.2 decimation 

Afin de simplifier la g6om6trie d'un maillage, il convient de d6cimer 

15 iterativement le maillage tout en conservant une bonne topologie, on choisit done 

une transformation elementaire : la fusion d' arete, qui fusionne deux sommets 
adjacents en un sonamet, supprime deux faces et place le nouveau sonraiet 
resultant de la fusion, suivant un principe decrit ult^rieurement. 

Cette decimation peut etre r^alisee altemativement ou successivement avec 

20 une optimisation des positions d'un ou de plusieurs sonunets suivant T application 

vis6e. On peut construire plusieurs niveaux de resolution d'un maillage en utilisant 

cette transfoffliatlon 616meniaire, el la giaiiulaiite niinimalc est obtenuc p ar une 

decimation puisque le maillage est alors defmi au sommet pres. On precede alors 
en alternant d6cimation locale et optimisation locale du sonunet resultant de la 

25 fusion. 

Inversement, on peut definir des niveaux de resolutions de maniere 
generale en decimant au prealable un ensemble de sonMnets puis en optimisant 
globalement les positions des sommets du maillage. 

La decimation successive permet de construire des maillages progressifs. 
30 Dans ce cas, il convient de choisir judicieusement Tordre de realisation des 
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transformations. Get ordre est defini par le biais d'une queue de priorite sur les 

transformations 61ementaires de maniere k efifectuer en premier les transformations 

affectant le moins possible le modele original. 

Cela amene h d6finir la notion de coflt pour chaque operation el^mentaire et 
5 la queue de priorite sera triee par ce demier. Ainsi initialement, toutes les 

transto rmaiioiis possibles sunt lislees, le cuul de chacune J'elles est calculg (saiiis 

effectuer de modifications svir le maillage), on realise la transformation de moindre 

coflt, et on recalcule les cofits des transformations elementaires modifiees par la 

transformation prec^dente. 
10 Bien sflr, d'autres techniques que la decimation peuvent etre utilisees pour 

obtenir le maillage approche, ou approxime, tel que des techniques de sous- 

6chantillonnage, de subdivisions adaptatives, etc... 

6.3. optimisation 

Le probleme de T optimisation consiste done k minimiser une fonctionnelle 
15 d' erreur definie par : 

E(M, M') = d(M, M') 
ou d(M, M*) caracterise une distance entre les deux maillages M et M% et garantit 
la fidelite aux donnees initiales. L'algorlthme d' optimisation est choisi notamment 
pour ses performances en terme de vitesse et/ou de pertinence de la solution finale. 
20 On explique ci-dessous le contexte theorique de la metrique basee volume 

utilis6e pour caracteriser Terreur d' approximation entre deux maillages, ainsi que 

Id formulation variationnellc pcrmettant Toptimisation de la quality de 

r approximation, selon 1' invention. 
6.3.1 metrique 

25 Un exemple de mesure de Terreur, ou distance, entre les maillages M et M' 

est illustree en figure 5, La mesure de 1' erreur est donnee par le volume V, 
compris entre les deux surfaces formees par les maillages M et M' . 
Ce volume est defini par T integrate de Lebesgue suivante : 

VCHIVT)- J Iv(M,M)(q)d'q 

R5 



• 
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5 


ox/*»/* T tn\ la frknrfinn \nr\\nstfr{n^ Hp 1'int^rieur du volume eto iin Doint 

gen6rique de . 

6.3,2 fonctionnelle 

Un maillage ^tant d^fini partiellement par les positions de ses sommets, on 
peut formuler le probldme de la minimisation du volume par le biais d'une 




fonctionnelle d energie a nunumser, cette roncuonneiie etant aeiinie par . 




E = Eeiror = M') = d(M.M'.Xi,...,X„)) 




L'algorithme d' optimisation mis en oeuvre consiste alors ^ determiner les 




positions optimales des sommets Xj , X2, ...Xn qui minimisent la distance 








O.J, J VCiriUUOn €l€TneTuCllr€ Cie vOl 




L» evolution eiemeniaire au maiiiage ivi esi ueiime par ic cnctuip uc 




vecteurs h^k . Une expression de la variation 616mentaire du volume est decrite en 




annexe 2 (§ 2). 
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6. 5. 4 mise en oeuvre du gradient adaptatif 




La technique d' optimisation utilis6e pour resoudre le probleme de la 




minimisation peut avantageusement reposer sur une approche iterative, tel que le 




gradient adaptatif. L' annexe 2 (§ 1) dterit en d6tail la mise en oeuvre de cette 


• 


technique du gradient adaptatif, ainsi que le calcul d'une variation 61ementaire de 
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volume permettant la determination du gradient & chaque iteration, et le 


25 


parametrage du maillage par un modele d'eiements tims tnangmaires. 

Un modele d' Elements finis triangulaires permettant de calculer le gradient 
en chaque sommet est d^crit en aimexe 2 (§ 3). 

L'algorithme du gradient pour un sommet Xj s*ecrit sous la forme : 

v., t 95d 



L' expression de la d6riv6e partielle de la distance pour un sommet Xi est : 
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d 5d 

— — = - JJ Tl(u,v)n(u,v)X.i(u,v)d<T(u,v). 

oXj u,v € Supp(A^ 

Le support de la fonction de forme est d^fini par les triangles voisins du 
sommet X|. L'integrale est calcul6e par le biais d'une discretisation sur les 

triangles (decrite par la suite) et une henristique a ete adoptee pour detenmnef 
5 Torientation relative des surfaces formees par M et M' (^galement decrite par la 

suite). 

A chaque iteration du gradient, on applique un deplacement sur le sommet 
Xi qui tend a minimiser le volume entre les surfaces. Le deplacement est la somme 
des contributions liees aux calculs d* orientation en plusieurs points de chaque 
10 triangle. 

Le pas Y- est adaptatif. Au depart il est initialis6 h (une fraction de la 

longueur de T arete fusionnee), puis il est multipli6 par k (avec k < 1), en presence 
d'une oscillation (c'est-^-dire lorsque deux vecteurs de ddplacement successifs 

sont de directions oppos6es : < — » ~dx7 ^ 

15 Ce pas Yj est born6 (borne inferieure Ymin^' ^^^^ conserver la ^ 

convergence de T algorithms Onay^G ]5, 2A,q-5], avec 0 < 5 < et II U^(k) 

II < 1/Xq (ou represente le Hessien de E). 

On considere la convergence atteinte lorsque le deplacement du sommet Xj 
d'une iteration k Tautre est inferieure a e (e etant une fraction de la dynamique 

20 locale des arStes du maillage). Le nombre d'iterations de Talgorithme peut 

^galement 8tre borne. 

7. description d'un mode de realisation detaille 

Le proced6 de T invention permettant de generer des approximations de 
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maillages a differents niveaux de resolution, un codage progressif de maillages par 
decimation et optimisation locale peut fetre mis en oeuvre. 

Dans cette optique, on definit une transformation 616mentaire de dteimation 
adapt6e au codage progressif, une methode d* optimisation locale et une queue de 
5 priorite sur les transformations afin de maximiser le ratio quantity 

d ' intofmalions/qualile visuelle a diaque ittiiaaun Jt ralgoiidunc . Oil JCuit pm la 
suite la mise en oeuvre de cet algorithme, et notamment le calcul du gradient d6fini 
sur un sommet du maillage, 
7. J transformation elementaire 
10 Le codage progressif permet de definir la notion d'6chelonnabilit6 sur un 

maillage. L'echelonnabilite maximale est obtenue lorsque la r6solution du maillage 
peut Stre choisi au sommet prds. Un maillage peut ainsi Stre code sous la forme 
d*un maillage de base (a une resolution visuellement acceptable) et d'une sequence 
de raffinements. 

15 Afin d'obtenir une telle structure de donnees, il convient de decimer 

iterativement le maillage, d'optimiser localement les positions des sommets de 
maniere a approximer localement au mieux le maillage initial, et d'enregistrer 
simultandment la sequence de decimation afin de coder le maillage dans le sens des 
raffinements. 

20 On choisit done une transformation elementaire reversible, par exemple la 

fusion d'arStes definie par HOPPE (document dej^ cite), qui est iUustree en figure 

1. La fusion d*ar£tc 10 consistc a fiisionncr les deux sonunets adjacents 11 et 12 

en un sommet 13, a supprimer les deux faces 14 et 15 et a optimiser la position du 
sommet 13 resultant de la fusion. 

25 Cette transformation est reversible (insertion 16 d'un sommet). 

Apres transfomiation, les adjacences doivent Stre mises ^ jour. Pour cela, 
on considere une structure de type faces/sommets et une adjacence definissant 
toute les relations entre les Elements. Les faces stockent des liens vers leurs faces 
voisines, et chaque sommet liste les faces qui le contiennent et les sonunets 

30 voisins. Les figures 12 et 13 illustrent la mise ^ jour des adjacences de faces apres 
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fusion d'une arete adjacente k deux faces (figure 12) ou h une face (figure 13). 

La decimation est ainsi naturellement prise en compte sur toutes les zones 
du maillage, sans qu'il soit necessaire de detecter les cas particuliers, telles que les 
arStes adjacentes a une seule face. 
5 ZiZ initialisation 

L' optimisation agit sur la position du gonunct resultant dc la fusion dc 
Tarete formee par les sommets et Xb- Afin de debuter Talgorithme 
d'optimisation avec une bonne condition initiale, on choisit de positionner le 
sommet le plus proche de roptimum probable. Pour cela, on introduit la notion 
10 d' arete vive : une arete est vive lorsque Tangle fonn6 par les normales aux deux 

faces adjacentes est sup6rieur k un seuU fixe parametrable. On denombre ensuite le 
nombre d'arStes vives autour des sommets de TarSte a fusionner, ainsi que cela est 
illustr^ en figure 14. Le sommet 141 est eissocie k aucune arSte vive, le sommet 
142 k deux ar6tes vives et le sonraiet a trois aretes vives. 
15 On deduit de ce denombrement deux cas pour T initialisation, ainsi que cela 

est illustre en figure 15 : 

si les nombres d' aretes vives autour de X3 et autour de X^ sont 

identiques, on effectue T initialisation au milieu du segment formant 
r arete a fusionner. Sur les zones planes de faible courbure, cela 
20 permet de conserver une bonne compacite sur les triangles voisins. 

Sur une arete vive reguliere (nombres d'aretes vives 6gaux a 2), 
cela permet de positionner le sommet proche de T arete vive qm sera 
conserve au cours de 1' optimisation ; 

si les nombres d'aretes vives autour de X^et autour X^sont 

25 differents, T initialisation est effectuee sur le sonunet presentant le 

plus grand nombre d'arStes vives. Dans les cas les plus courants, 
roptimum est atteint a partir de cette position initiale. 
Dans Texemple de la figure 15, dans lequel le maillage source correspond 
a un parallelepipede 151, on constate que cette heuristique place le sommet : 
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sur le coin du parall616pipede lorsque Tarete h fusionner forme un 
coin 152 ; 

et sur 1' arete vive reguliere du parall61^pip&de lorsque I'ardte 
commence sur la region plane et se termine sur Tarfete vive 155. 
Les situations oil le nombre d' aretes vives est le meme autour des deux 



sommets est illustrd en 152 et 154. 
7.3 optimisation 

On rappelle que T expression de la fonctionnelle consid^ree : 

E = Ee™r = d(M,M') 

Cette fonctionnelle est minimis^e en optimisant la position du sonmiet 
par un algorithme de gradient : 

asd' 



Les determinations du pas et des critferes d'arrSt sont discutees au 



§ 6.3.4. 

15 7.4 discretisation 



Le calcul du gradient necessite done une discretisation de I'int^grale : 
-;j Ti(u,v)n(u,v)>.j(u,v)da(u,v) 

U.V e Supp(X^ 



Sur le support de la fonction de forme d6fmie sur Xi. 
On est done amen6 a dchantillonner la surface des triangles, avant 
20 d'evaluer en chaque point X(x,y,z) la valeur de la fonction de forme X,.(x, y, z) et 

r orientation des surfaces 'n(u,v). 

Le gradient s'exprime done sous la forme discrete suivante : 

asd 

3.r=-^ I L r|(i, J,)S(i„jjA,,(i„jx)da(i„jJ 
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avec : - S ^ensemble des triangles voisins du sommet . 

le nombre de points echantilloimes dans la direction de u ; 
le nombie de points echantillonnes dans la direction de v. 

Po\ar un triangle donne, rechantUlonnage est effectue dans le plan passant 
5 pax 3ommctG, le nombre de points rt-^nt proportionnrl a snrfanr. dn triangle, 

ainsi que cela est illustre en figure 16. 

Le repere d'6chantillonnage est forme par 1' arete 161 de plus grande 
longueur et la hauteur 162 du triangle courant. Soit n le nombre minimum de 
points par triangles fixes a r avance, S^^ I'^ire minimale des triangles voisins de 
10 Xi et h le pas d'echantillonnage. On deduit h par la formule suivante : 

La fonction de forme est evaluee par un rapport de surfaces tel que decrit 
en annexe 2. 

7.5 orientation des surfaces 
15 L' orientation relative de deux surfaces definie en trois dimensions est 

delicate a d6finir avec exactitude. On utiUse avantageusement une heuristique pour 
calculer r orientation des deux surfaces M et M' (c'est-a-dire le terme r|(u, v) de 
Tequation du gradient present6 ci-dessus) : 

r|(u,v) = -<ni^,nM>.<aM', XjX^ > 

20 avec : Xj le point echantillonne sur M' ; 

X le point d' intersection de la droite passant par X. et de direction fij^ 

avec le maillage source M ; 

la noraiale au maillage source M au point Xj^ ; 

n^. la noraiale au maillage source M' au point Xj. 

25 ou <,> repr6sente le produit scalaire. 

La figure 17 illustre ce calcul de Torientation en deux dimensions. En trois 



• 
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dimensions, il suffit de remplacer les segments par des triangles. Cette heuristique 
limite fortement la foraiation de plis sur la surface, puisque le terme < Bm.3m> 

traduit 1' opposition des normales lorsqu'elles sont de sens opposes. 
7,5 priorites sur les transformations 
5 Chaque transformation 61^mentaire d^cime le maillage approximant M' . La 

quality de T approximation se degrade done au cours de la decmiation, ou reste au 
mieux invariante lorsque le volume reste inchang6, et ce malgr6 Toptimisation des 
points. 

Par exemple, si Ton considere une sphere, on sait qu'elle necessite en 

10 theorie line infinite de triangles pour Stre parfaitement mod61isee. Des lors que Ton 

restreint le nombre de triangles, on degrade la qualite de son approximation. En 
revanche, sur des regions planes d'un objet, on peut d^cimer jusqu*& Tobtention 
d'un maillage optimal (ainsi, sur un parall616pipfede, T optimum est atteint lorsque 
chaque face est decrite par deux triangles). 

15 On souhaite coder progressivement les maillages, et done retirer 

it6rativement les sommets par une op6ration de fusion d' aretes. Aprfes 
r optimisation des positions, un autre processus est done d6velopp6, consistant & 
definir une priority sur les transformations eldmentalres, de maniere ^ eff ectuer en 
premier les transformations affectant le moins possible le mod&le. 

20 Pour cela, on definit la notion de coflt pour chaque operation 61ementaire, 

et une queue de priorite triee selon ce coftt. 

Aittsi, iriitialement, toures les transformations possibles sunt list^es puis le 

coiit de chacune d'elles est calculi (sans effectuer les modifications 
correspondantes sur le maillage). On realise la transformation de moindre coOt, 

25 puis on recalcule les coflts des transformations elementaires modifiees par la 

transformation pr€c6dente. On rappelle que la fusion d'arSte supprime deux faces 
soil en consequence, les transformations concemant ces faces sont retirees de la 
queue de priorite, ce qui cree une ou deux nouvelles transformations, pour 
lesqueUes on calcule 6galement le coflt en Anergic, deduit de la technique decrite en 

30 annexe 2, et qui repr^sente la variation de volume entre le maillage avant et aprfes 



24 



transformation. 

Le calcul direct du volume entre deux surfaces triangul6es 6tant complexe k 
calculer, on fait appel a une heuristique permettant de simplifier ce probleme. 
Soient : 

5 - Tj la transformation fusionnante de sommet Xj et Xj de M' ; 

Xj^le sormnet de M' resultant de la fusion (dont on a optimise la 
position) ; 

F(X/) les faces de M' avoisinant le sommet X^^ apres 
transformation ; 

10 - I'ensemble des sonmiets de M appartenant aux faces ayant 6te 

intersectees lors du calcul de 1' orientation des surfaces pendant 
r optimisation. 

Le coflt d'une transformation 61ementaire Tj s'exprime alors sous la 

forme : 

15 C(Ti(Xi,Xp) = maxd2(VM,F(Xi^)) 

Ce coflt correspond done h la distance maximale des sommets du maillage 
original (sous les faces avoisinant le sommet Xj) aux faces avoisinant Xj apres 
transformation. La figure 18 illustre, en deux dimensions, I'expression du coflt 
dans le cas d*une decimation sur une discontinuity. Comme indique 

20 precedemment, en trois dimensions les segments sont remplaces par des triangles. 

La transformation 181 consiste a fusionner les sommets Xj et Xj en un sonunet 
Xj'. Les sommets 182i a I824 correspondent k Tensemble V^. 

On peut constater que sur un plan maille, le coflt d'une transformation 
serait nul, les zones comportant une information redondante etant d6cimees en 

25 priorite. 

On presente maintenant des exemples de synoptique pour la mise en 

oeuvre des diff6rents aspects discut6s precedemment. 

7,8 synoptiques 



25 





7. 8. 1 approximation de maillages 

La figure 19 resume Talgorithme d' approximation d'un maillage M 
comnortant n faces oar un maillage M' comportant m faces (avec bien s(ir m < n). 

On conmience tout d'abord par copier le maillage source M, ou maillage 
original, dans la variable M', reprdsentant le maillage simplifi6 que Ton cherche h 




determiner (191). Sur ce maillage M', on effectue ensuite une recherche de toutes 




les transformations dlementaires r^alisables (192), puis un calcul (193) des codts 




6nerg6tiques de ces transformations aprfes optimisation. 




On s61ectionne et on realise (194) la transformation de moindre coflt parmi 
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Tensemble des transfomiations possibles. D s'agit de Fetape de decimation. 




Puis, on met i jour (195) les coflts 6nerg6tiques des transformations, dans 




le voisinage affectd par la transformation effectuee h I'^tape 194. Enfin, on ajoute 




(196) les nouvelles transformations induites par la transformation effectuee 194, et 




on calcule les coflts correspondants. 
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Les trois 6tapes 194, 195 et 196 sont r^it^rees (197) n - m fois, jusqu'k ce 




qu'au moins un des criteres d' arret soit atteint. 




7. 8, 2 optimisation en position 




La figure 20 illustre T optimisation du sommet Xi, aprfes fusion de TarSte 


m 


(Xi, Xj). 


20 


On definit au prealable les paramdtres d' initialisation suivants. 




partiellement en fonction de la dynamique de la geometric locale autour du sommet 
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Xi : 

Epsilon est un critere d'arret de la convergence, fraction de la 
longueur moyenne des aretes autour du sommet Xj (typiquement 

0,001) ; 

Maxlter est la borne superieure du nombre d' iterations, fix^e 
arbitrairement (typiquement 200 iterations) ; 

Steplnit, StepDec et StepMin sont respectivement la valeur 
d'initialisation, le taux de d6croissance et la bome inf6rieure du 
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pas. Steplnit et StepMin sont des fractions de la longueur de T arete 
a fusionner (typiquement 0,1 et 0,01), StepDec est fixe 
arbitrairement (typiquement 0,95) ; 

StepSampling est le pas d*echantillonnage sur les triangles, utilise 
pour le calcul discret de I'int^grale formant le gradient a chaque 
iteration. H est ddduit des aires des faces voisines au sommet A. 
NbSharpEdges(X) est le nombre d' aretes vives autoxir du sommet 
X. 

L'algorithme debute par une etape 201 d' initialisation en position, qui 
consiste a effectuer les operations suivantes : 

Si (NbSharpEdge(Xi) = NbSharpEdge(Xj)) 

alors 

initialiser au milieu de 1' arete a fusionner, 

sinon 

initialiser sur le sommet du plus grand nombre d'aretes vives. 
On realise ensuite le calcul 202 du gradient, avec les operations de 
discretisation sur les triangles et d' orientation des surfaces. Puis, on effectue le 
deplacement 203. Le sorrunet se deplace de la valeur du vecteur gradient multiplie 
par le pas courant. 

On contr61e ensuite revolution 204 du pas. Si une oscillation est detectee 
avec le vecteur gradient precedent, le pas est multiplie par StepDec. Si le pas est 
inf^rieur a StepMin, il est reinitialise a StepMin. 

Enfin, on effectue un test de convergence 205. Si la distance entre la 
position courante et la position precedente est inferieure a Epsilon, ou si le nombre 
d'iterations Maxlter est atteint, Talgorithme est arr8te (206). Sinon, on repete 
(207) les etapes 202 k 205. 
7.9 resultats 

On montre maintenant, k partir d'un exemple simple (un peirallelepipede), 
le comportement du proced^ d' approximation de maillages selon Tinvention. 
L'interet du parallelepipede est de permettre de verifier la qualite de 
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r approximation sur les discontinuites (arStes vives r^gulieres et coins) et sxir les 
plans. 

Sur un tel objet, 1' optimum est atteint lorsque Ton obtient chaque face 
maill6e par deux triangles et que la forme globale de T objet est parfaitement 
5 conserv^e. 

La flgUA ^ fe 21 illustre le cas d'une ddcimalion hu t un coin 211 dtp= 
parall616pipfede. La solution initiale a 6t6 plac6e sur le coin 211, puisque le nombre 
d*ar6tes vives est maximal en ce sommet (trois autour du coin 211 et deux autour 
de TarSte viva reguliere de Tobjet (sommet 212)). L'algorithme d' optimisation 
10 oscille, puis converge autour de cette position initiale 213, puisque le pas de 

gradient a ete initialise k une valeur non nuUe, puis a diminue dfes que Toscillation 
a 6t6 d6tect6e. 

La figure 22 presente une decimation sur une arSte vive reguliere. La 
solution initiale a 6t6 plac^e au milieu de TarSte 223, puisque le nombre d' aretes 

15 vives est identique (deux autour de chaque sommet 221 et 222). Plusieurs 

solutions situ^es sur TarSte vive reguliere de Tobjet satisfont Tinvariance de 
volume, et Talgorithme converge apres oscillation autour de TarSte, sur une 
position d^pendante de la solution initiale. Le milieu 224 de TarSte constitue done 
une bonne initialisation dans ce cas. 

20 La figure 23 illustre une autre situation de decimation sur une arete vive 

r^gulifere, dans laquelle la solution initiale a 6t6 placfe sur TarSte vive 234, puisque 

les nombrea d'ordtcs vivos sont diffdrcnts (deux autour du sommet 233 situ6 sur 

r arete, et aucune autour du sonmiet 232 situe sur la face plane 235 de Tobjet). 
Plusieurs solutions situ&s sur T arete vive reguliere de T objet satisfont Tinvariance 

25 de volume, et l'algorithme converge aprds oscillation autour de TarSte 234, sur 

ime position 236 dependant de la solution initiale. 

La figure 24 pr€sente le cas d'une decimation sur un plan 243. La solution 
initiale a 6t6 placee au milieu de TarSte 245 formee par les sommets 241 et 242, 
puisque les nombres d'arStes vives sont identiques (aucune autour de chaque 

30 sonunet 24 let 242). Plusieurs sommets situds sur le plan de T objet satisfont 
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rinvariance de volume. L'algorithme converge, apres oscillation, de part et d* autre 
du plan, sur une position dependant de la solution initiale. Le milieu 244 de TarSte 
constitue done, a nouveau une bonne initialisation dans cette situation. 

Le procede decime le maillage initial en utilisant une queue de priority. 
5 La figxire 25 montre les douze iterations successives de l'algorithme, 

numeio Lees (b) a ( m), a pailii du maillage oilgiiial uum^ioLtS (a). 

On constate que les arStes entrainant une invariance de volume sont 
decimees en priority. L*optimum est atteint pour les sommets demandes sur cet 
objet, puisque chaque face est alcrs maillee par seulement deux triangles. 
10 On remarque que les aretes rejoignant deux coins de Tobjet ne sont pas 

fusionndes, puisque cette operation aurait entralne une forte degradation de la 
quality de I'approximation, les couts et transformations elementaires 
correspondantes dtant superieurs aux couts d*une fusion d'arStes sur un plan ou le 
long d'une arete vive r6guliere. 
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ANNEXE 1 



^ Codage progressif de maillages 

Huguta Huppe adevelupp^ uiie iiieLliude fuucLiunnant par decimation et optimisation d6s 
points. La decimation est realis^e pax une transformation elementaire reversible et facile- 
ment codable : la fusion d'arete qui rassemble deux sommets adjacents en un seul sommet, 
supprime deux faces et optimise ensuite la position du sommet resultant de la fusion (figure 

1). 

L'auteur a utilise une approche variationnelle permettant la minimisation du proces- 
sus, la fonctionnelle d'^nergie associee est d^finie par un terme d'erreur et un terme de 
YQ regularisation : 

E = E^rror + Ercgul 

Ce terme d'erreur caracterise la distance entre les deux maillages M et M definie lo- 
calement, et garantit la fidelite aux donnees initiales. Le terme de regularisation permet de 
respecter la topologie d'un maillage et d'assurer I'unicite de la solution. Afin d'initialiser 
le probleme, une fois la face a supprimer determinee, trois minimisations (partant de trois 
positions initiales differentes choisies de maniere a initialiser correctement le processus d'op- 
timisation) sont realis^es (figure 2) afin de retenir la solution correspondant au minimum 
25 d'energie de la fonctionnelle. 

Chaque transformation Ti supprimant un sommet de la triangulation et modifiant le 
maillage Mi en Mf+i , engendre une augmentation de I'energie globale qui peut etre calculee. 
Le cout en energie resultant de cette decimation est le suivant : 



{ 



AExi = £?jifi+i - EMi 
AExi > 0 



Au depart toutes les transformations possibles sur le maillage sont listees, puis le cout 
2Q energ^tique pour chacune d 'entre elles est calcule. Ces transformations sont rangees dans une 

queue de priority triee suivant leur cout estimd. Au cours de Talgorithme, la transformation 
de cout minimum est effectu^e, puis les couts des transformations dans le voisinage afFecte 

pai cbtle inodiflcatlou soul mis ajuui (figuie 8). Les liajisfoimaliunb iiuuvdleinent cr^fas sunt 

ajoutees et leur cout estime est calcuie (figure 4). L'algorithme agit done de maniere iterative 
et decime le maillage M d'un sommet et de de\ix faces a chaque iteration. Cette approche 
de decimation et d 'optimisation est bien adaptee au contexte de codage progressif, puisqu'il 
suffit d'enregistrer la sequence des decimations dans le sens inverse (celui des raffinements) 
pour coder un maillage sous une forme progressive. Ainsi un maillage M peut etre represente 
25 par un maillage a la resolution Mq et par une sequence de raffinements {ro,ri,r2, - . . , r„}. 

L'erreur utilisee pour caract^ser la distance entre les deux maillages est la distance 
euclidienne (ou quadratique) entre le point Xi a optimiser du maillage M et le maillage M, 
qu'il convient de definir. Les points situes a une distance D de Xi dans forment une 
sphere de rayon D : 

d2 (Xi, M) = min da {Xi, Y) = min HAT,- -Y\\^ (2) 



30 



30 



ou 

5 ll-Y||2 = (f;x?)^/^x€K^ (3) 

Le maillage M etant une surface triangulee, cela engendre une mesure de la distance entre 
un point Xi et un ensemble de triangles. On remaxque d'apres la definition de la distance 
dans que le point de M situe a la distance minimale de Xi n'est pas toujours unique. Le 
choix effectif de ce point depend de Tordre d'insertion des points et de la nature du test lors 
de la recherche du minimum (< cm <). La metrique appliquee sur les maillages devient 
la metrique puisque Ton raisonne avec une distance point-surface. 

10 
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20 
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ANNEXE 2 



Optimisation 



5 La technique d*optimisation choisie pour resoudre notre probleme de minimisation est 

une methode iterative, appelee gradient adaptatif. 

1 Gradient adaptatif 

Notre objectif consiste done a minimiser la fonctionnelle E decrite precedemment. Comme 
dans tous les cas, la stationnarite de E est une condition necessaire d'optimalite, effective- 
ment, pratiquement toutes les m^thodes d'optimisation sans contraintes dans consistent 
10 a rechercher un vecteur X stationnaire (VJS;(A') = 0). Ce probleme est equivalent a la 

resolution du systeme d 'equations non lin^aires : 

^(x) = o,vi = i,...,;v (1) 

On peut chercher a resoudre directement ce systeme, ce qui conduit a la methode de Newton. 
Cependant, cette methode peut ne pas converger, si le point de depart des iterations est trop 
15 eloigne de A'. D'autre part, elle suppose la fonction deux fois continument difFerentiable, et 

elle necessite le calcul des derivees secondes en chaque point. 

C'est pourquoi les methodes les plus couramment utilisees precedent difFeremment : il s'agit 
de procedures iteratives ou Ton engendre une suite de vecteurs A'°, A^ ... convergeant vers 
un optimum local de E. 

Ces methodes de gradient constituent une famille de methodes qui precedent de la fa^on 
suivante : On part d'un point et on calcule le gradient V£;(A'°) en A"° de la fonction E 
a minimiser. Comme VEiX°) indique la direction de la plus grande augmentation de E, on 
se deplace d^une quantite 70 dans la direction opposee au gradient, et on definit le point : 

yl - vo , ^ (2) 
^ '^'""WVEiXO II 

La procedure est repetee et engendre les points A'°,XS A:",... suivant la relation : 

Dans cette famille, il convient de signaler les methodes de gradient a pas determine 
dans lesquelles on choisit a priori les vaJeurs de deplacements 7^. L'inconvenient de cette 
procedure est que la convergence peut etre tres lente. Nous avons done choisi d'utiliser la 
methode de gradient a pas adaptatif, qui consiste a diminuer (resp. augmenter) le pas 7fe 
lorsque I'erreur augmente (resp. diminue), tout en conservant la contrainte de controle du 
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pas, soit : jk € [jmin^Jmax]- Finalement, comme critere d'arret, nous nous sommes fixes que 
la difference entre deux erreurs successives devait etre inferieure a un £ donne. 

rn "-'^ rlr Iff minimiRation. les positions des sommets du maillage vent ^voluer de 
maniere a minimiser le volume entre les deux surfaces, cette evolution necessiie dts liaiculex 
le gradient a Titeration k (le terme VE{X^) de Tequation 

2 Variation elementaire de volume 

Le processus d 'optimisation est done vu comme une evolution du maillage M de maniere 
a minimiser V'(A/,M). L'evolution elementaire du maillage est definie par le champ de 
vecteurs SKl. En realisant une parametrisation de la surface par u et v, le vecteur 6M 
pent s'exprimer sous la forme Ji&(u, v). Si on approxime localement la surface par son plan 
tangent la variation elementaire de volume est ajors le parallepipede genere par revolution 
de I'element de surface dudv dans la direction SM{u,v) (figure 6). 

On peut ainsi exprimer la variation elementaire de volume induite par I'element d'aire 
dudv : 



6V{M,M,6M{u,v)) = 



dM{u,v) dM{u,v) 
du dv 



• SM{UyV) 



• dudv 



La variation est negative lorsque le vecteur 6M est oriente vers I'interieur du volume 
d'erreur. la figure 7 presente le cas particulier ou revolution du maillage se fait dans la 
direction de la normale. Prenons maintenant en consideration le fait que : 



du dv \ 



• dudv = n(u, v)da{u, v) 



avec n{u,v) la normale unitaire en M{u,v) et da^u^v) = ||^ A^W'dudv I'element d'aire. 
On exprime la variation de diotance entre M et M : 



6d{M,M,SM) = f f \n(u,v)6K<[{u,v)^da{u,v 



) 



Cette expression n'est pas facilement utilisable a cause de la valeur absolue. Nous intro- 
sons maintenant la fonction : 7j(u,«) = r]{M, M,fHu,v)) = 1 si la normale est orient^e 
s le volume interieur et -1 sinon, nous obtenons alors 

|n(u,t;)<5j^(u,u)| = »j(ti,i;)n(u,u)5M(u,v) 
6diM,M,S^) = - f f v{u,v)n{u,v)6M{u,v)da{u,v) 

J JutV 



d'ou 
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3 Modele d'elements finis triangulaires 

Afin de cadculer le gradient en chaque sommet du maillage a Titdration k de ralgorithme 
d'optimisation, on parametre le maillage par un module d'^l^ments finis triangulaires : 



N 



(4) 



-Yt est un sommet du maillage defini dans ^ et Xi{u, v) est une fonction de forme adaptee a 
un modele d'elements finis triangulaires. On definit une fonction de forme sur le maillage par 
le biais de coordonnees barycentriques tel que cela est decrit dans le document "la methode 
des elements finis, formulation de base et problemes lineaires", par O.C. Zienkiewicz et R.L. 
Taylor (afnor technique, 4ieme edition, 1991). Pour cela on numerote les sommets du triangle 
(figure 8), et on definit un systeme de coordonnees Ai, A2 et A3 liees par les relations : 

X = Ai • xi H- A2 * 3:2 + A3 • X3 

y = Ai • 2/1 + A2 • y2 + A3 • j/3 
1 = Ai -h A2 + A3 

Nous allons maintenant determiner les valeurs des fonctions de base : au sommet 1, 
Ai = 1, A2 = 0 et A3 = 0. Les lignes de niveau de Ai sont des lignes equidistantes paralleles 
au cote 2 — 3 le long duquel Ai = 0. On peut exprimer la coordonnee Ai en un point P par 
un rapport d'aires : 

A, = (5) 



at re (123) 



On obtient : 



^ _ ai -h bi • X -h Ci 


•y 




2A 








^ ^2 + 62 * X -h C2 


•y 




2A 










•y 




2A 










1 


Xi 


yi 


azre(123) = | ' 


1 




yt 




1 


X3 


yz 
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Avec : 



ai = X2 • ys - • y2 
61 = y2 - ys 

Ci = X3 — X2 

On calcule les valeurs 02, 62) C2, 03, 63, et C3 par permutation circulaire des indices 1, 2 et 3. 
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4 Evolution de surface 

L'^volution du maillage peut ainsi s'ecrire 



5d(M , M, SXi , SXi, .... SX„) = -J2 ^)"(". v)SXiXi(u, v)d(T(u, v) 

on en deduit la d6rivee partielle de la distance pour un sommet Xi : 
dV f f 

J Ju.veSuppiXi) 

L'aJgorithme de gradient pour un sommet Xi s'ecrit sous la forme : 



'^''^'dxl 

ou 7,- represente un pas adaptatif avec 7^ 6 [7min,7o]. Au depart le pas est initialise a 70, 
puis i^l^est multiplie par k (avec A: < 1 ) en presence d'une oscillation, c'est a dire lorsque 
22 ^ ^eXi » Wc7 X 0 Le comportement de I'algorithme d 'optimisation peut-etre illustre en 2d 

(figures 10 et 11), revolution du sommet Xi du modele M est liee a la somme des vecteurs 
d'orientation du modele par rapport a la courbe originale M, pond^res par la fonction de 
forme qui vaut 1 en Xi et decroit en 0 le long des segments voisins de Xi. En 3d le principe 
est le meme avec une fonction 'de forme definie sur les triangles voisins du sommet Xi. 
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REVENDICATIONS 

1 . Proc6d^ de codage d'un maillage source (M) representatif d'un objet en 
trois dimensions, d61ivrant un maillage simplifie (M*) correspondant audit maillage 

5 source (M), lesdits maillages 6tant d6fmis par un ensemble de sommets, d' aretes 

et/ou de faces, 

caract6ris6 en qu'il met en oeuvre une 6tape de minimisation d'un volume compris 
entre ledit maillage source (M) et ledit maillage simplifie (M'). 

2. FToc6d6 de codage d'un maillage source selon la revendication 1, 
caracterise en ce que, chacun desdits maillages 6tant d6fini par la position de 
chacun de ses sonraiets, ladite 6tape de minimisation assure la determination des 
positions des sommets (X j, Xj,... X^) dudit maillage simplifie (M') minimisant le 
volume V(M, M') entre ledit maillage source et ledit maillage simplifie. 

3 . Pxoc6d6 de codage d'un maillage source selon I'une quelconque des 
revendications 1 et 2, caract6ris6 en ce que ladite etape de minimisation met en 
oeuvre un processus iteratif optimisant progressivement les positions des sommets 
dudit maillage simplifie (M'). 

4 . Proced6 de codage d'un maillage source selon la revendication 3, 
caracterise en ce que ledit processus iteratif est interrompu lorsqu'au moins un des 
criteres d' arret suivants est atteint : 

^ un nombre maximum d'it6rations ; 

une difference entre deux vecteurs d^placement successifs des 
positions des sommets inferieure k un seuil predetermine (Epsilon). 

5 . Procede de codage d'un maillage source selon Tune quelconque des 
25 revendications 2 a 4, caracterise en ce que ladite etape de minimisation met en 

oeuvre une methode de gradient adaptatif . 

6. Procede de codage d'un maillage source selon la revendication 5, 
caracterise en ce que ladite methode de gradient adaptatif repose sur les operations 
sviivantes : 
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selection d'un vecteur Xp de R ^" (n > 1) dudit maillage simplifie, 

et calcul du gradient VE (Xp) en Xp de la fonction a niinimiser 
E = d(M, M'(Xi,...,X„)) ; 

d6termination de la position X^,* de Xp dudit maillage, salon la 


5 


relation definie a I'iteration k+1 par : 




VE(X„'') 

X ''■"-X "-Y.x „ ^ „ 




^"^^ ||VE(X/)|| 




k variant de 0 a n-1 et y^^ etant le pas de ladite relation. 




7 . Proced^ de codage d'un maillage source selon la revendication 6, 




caract6ris6 en ce que ledit pas est variable, et en ce qu'il varie en fonction de 


10 


r oscillation de deux deplacements successifs des sommets et/ou en fonction des 




variations de I'energie. 




8 . Precede de codage d'un maillage source selon la revendication 7, 




caract6rise en ce que ledit pas est maintenu entre deux bomes et . 




9 . Proc6d6 de codage d'un maillage source selon Tune quelconque des 


15 


revendications 2^8, caract6rise en ce que, a chaque iteration, on determine une 




variation elementaire dudit volume correspondant k un champ de vecteurs 5^ , 


i 

20 


et en ce que, la surface etant parametree par u et v, de fagon qu'un vecteur Sf^ 
s'exprime sous la forme SS/I* (u, v), ladite variation el6mentaire est assimilee au 
parallel6pipede g6nere par revolution de T element de surface dudv dans la 
direction SSV (u, v). 

10. Precede de codage d'un maillage source selon la revendication 9, 
caracteris6 en ce que la variation elementaire entre M et M' s'6crit : 
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5d(M,M",5S*) = - / J ii(u.v)n(u,v)8Sl'(u.v)da(u,v) 

avec: r|(u,v)=r|(M,M,n(u,v)) = 1 si la normale unitaire n(u,v) est 

orient^e vers le volume intdrieur, et -1 sinon. 

1 1 , Prnc.M^. de codage d*nn maillagft source selon Tune quelconque des 
5 revendications 1 k 10, caracterise en ce que ledit maillage simplifie est parametre k 

Taide d'un modele d' Elements finis. 

12. Proced6 de codage d'un maillage source selon la revendication 11, 
caract6ris6 en ce que lesdits elements finis sont obtenus & Faide d*un interpolateur 
affine. 

10 13. Procede de codage d*un maillage source selon Tune quelconque des 

revendications 1 1 et 12, caract6ris6 en ce que ledit maillage simplifie s'6crit : 

N 

M(u,v)=5:Xi>.i(u.v) 

i = l 

avec A,.(u,v) une fonction de forme adaptee h un modele d'el6ments finis 
triangulaires, 

15 et Xj sommet dudit maillage, defini dans R^, 

et en ce qu'on d^finit une fonction de forme sur le maillage par le biais de 
coordonnees barycentriques. 
TT: Procede de codage d'Urt maillage source seluii la leveudication 13, 

caracterise en ce que la variation de la position d*un sommet Xj ^ la k'^'"^ iteration 
20 s'dcrit: 

ic+i k 5 5d 



ax, 



avec la deriv6e partielle de la distance pour un sommet Xj : 
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— — = - J J Ti(u,v)n(u,v)^i(u,v)da(u,v). 

o^i u,v eSupp(Aj 

15. Precede de codage d'un maillage source selon Tune quelconque des 
revendications 1 h 14, caracterise en ce qu'il met en oeuvre un codage progressif 
dudit maillage simplilie, par decimation d'lin sommet et optimisation locale ou par 

5 decimation d'un ensemble sommets et optimisation globale. 

16. Precede de codage d'un maillage source selon la revendication 15, 
caracteris^ en ce que ledit maillage simplifie est represent^ par un maillage de base 
et une sequence de raffinements dudit maiUage de base. 

17. Precede de codage d'un maillage source selon Tune quelconque des 
10 revendications 15 et 16 et selon la revendication 6, caracteris6 en ce qu'il 

comprend une etape d'initialisation du veceteur Xp, comprenant une etape de 

denombrement des arStes vives autour de deux sommets a fusionner, apr^s ladite 

decimation, 

et en ce que : 

15 - si les nombres d'arStes vives sont les memes autour des deux 

sommets, on place le soiranet resultant de la fusion au milieu du 
segment reliant lesdits sommets ; 

si les nombres d' aretes vives sont differents, on place le sommet 
resultant de la fusion sur le sommet presentant le plus grand 
"20 nombre d' aretes vives. 

18. Precede de codage d'un maillage source selon Tune quelconque des 

revendications 14 a 17, caracterise en ce que le calcul du gradient a k>^"^^ iteration 

comprend les etapes suivantes : 

discretisation de 1' expression de la derivee partielle de la distance 
25 pour chaque sommet Xj, sous la forme : 

- o J N M 

5 5d " ' 

■T^=-E E L'n(ixJx)n(i.JT)>^i(itJx)dCT(i^Jx) 

TeSi=lj=l 
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avec : - S Tensemble des triangles voisins du sommet Xj . 

le nombre de points 6chantillonn6s dans la 

direction de u ; 

le nombre de points 6chantillonn6s dans la 

direction de v ; 



orientation des surfaces, par identification de la plus proche 
intersection avec le maiUage source M. 

19. Proced6 de codage d*un maillage source selon Tune quelconque des 
revendications 1 a 18, caracteris^ en ce que 1' orientation relative des surfaces audit 

10 maillage source et dudit maillage simplifi6 est calculi selon T Equation : 

n(u,v) = -<nM.,nM>.<n^, XjXM> 

avec : Xj le point dchantillonnfi sur M' ; 

Xj^ le point d* intersection de la droite passant par Xj et de direction n 

avec le maillage source M ; 
15 la normale au maillage source M au point Xj ; 

n^. la normale au maillage source M' au point Xj. 

20. Proced6 de codage d'un maillage source selon Tune quelconque des 
revendications 1 a 19, caracterise en ce quMl comprend une 6tape de limitation de 
la degradation due h une transformation 616mentaire, mettant en oeuvre une queue 

"20 de priorite sur les transtormations eiementaires. 

21. Proced6 de codage d'un maillage source selon la revendication 20, 
caracterise en ce que ladite etape de limitation de la degradation due a une 
transformation 61ementaire, c'est-^-dire une fusion d'arSte, definie par deux 
sonunets, comprend les 6tapes de : 

25 - calcul d'un coflt pour chaque transformation eiementaire possible ; 

realisation de la transfomiation elementaire de moindre coflt ; 
recalcul des cofits des transformations 616mentaires modifiees par la 
transformation elementaire precedente ; 
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ajout des nouvelles transformations 61ementaires creees et calcul 

des coiits correspondants. 
22, Procede de codage d'un maillage source selon la revendication 21, 
caracterise en ce que le coflt d'une transformation elementaire (Tj) s'exprime par : 

C(T,(Xi,Xp) = maxd2(V^,F(X[)) 
avec : Tj transformation fusionnant deux sommets Xj et Xj dudit maillage 

simplifieM' ; 

X/ le sommet dudit maillage simplifie M' resultant de ladite 
transformation ; 

F(X|0 les faces dudit maillage simplifie M* avoisinant le sommet X-^ apres 

ladite transformation ; 
Vm ensemble des sommets dudit maillage source M appartenant aux 

faces ayant 6t6 intersectees lors du calcul de T orientation des 

surfaces pendant ladite minimisation. 
23. Application du proc6de de codage d*un maillage source selon Tune 
quelconque des revendications 1 k 22 i au moins un des domaines appartenant au 
groupe comprenant : 

la realite virtuelle ; 

la simulation scientifique ; 

la modelisation. 



1/11 




2/n 




Fig. 7 




Fig. 10 



Fig. 1 1 



4/n 




Fig. 1 3 



142 — — ><(\ 


A 13 


141 


^ — 


V3 





Fig. 1 4 




Fig. 1 5 



6/11 




Orientation relative 
des surfaces 



Fig. 1 7 



7/n 

Xi 





> { 




• 


i 




< 


> • 



• Sommets d fusionner 



^^^^"'^ max\ 


Xi 




c? o • • 



© Vm 

Fig, 1 8 



8/11 



191 



192 



193 
197 



(n-m) 



Fig. 19 



194 

195 
196 



207 



Copie : M — > M' 
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et calcul des coOts correspondants 
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Calcul du gradient 
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Fig. 20 



Test de convergence 
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Fig, 22 
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(a) Module original 



(b) Iteration 1 



(c) Iteration 2 




(d) Iteration 3 (e) Iteration 4 (f) Iteration 5 
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Fig. 25 
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